arXiv: 1509.01914v2 [cond-matsupr-con] 20 Nov 2015 


Josephson current in a normal-metal nanowire coupled to 
superconductor/ferromagnet/superconductor junction 


Hiromi Ebisu^, Bo Lu^, Katsuhisa Taguchi^’^, Alexander A. Golubov^’^, Yukio Tanaka^’^ 

^ Department of Applied Physies, Nagoya University, Nagoya fGf-SGOS, Japan 
^ Department of Physies, Hong Kong University of Seienee and Teehnology, Clear Water Bay, Hong Kong, China 
^ Faeulty of Seienee and Teehnology and MESA-h Institute for Nanoteehnology, 

University of Twente, 7500 AE Ensehede, The Netherlands 
^ Moseow Institute of Physies and Teehnology, Dolgoprudny, Moseow 141700, Russia 

(Dated: November 23, 2015) 

We consider superconducting nanowire proximity coupled to superconductor / ferromagnet / 
superconductor junction, where the magnetization penetrates into superconducting segment in 
nanowire decaying as ~ exp[—^], where n is the site index and the ^ is the decay length. We 
tune chemical potential and spin-orbit coupling so that topological superconducting regime hosting 
Majorana fermion is realized for long We find that when ^ becomes shorter, zero energy state at 
the interface between superconductor and ferromagnet splits into two states at nonzero energy. Ac¬ 
cordingly, the behavior of Josephson current is drastically changed due to this “zero mode-non-zero 
mode crossover”. By tuning the model parameters, we find an almost second-harmonic current- 
phase relation, sin 2(p, where (p is the phase difference of the junction. Based on the analysis of 
Andreev bound state (ABS), we clarify that the current-phase relation is determined by coupling 
of the states within the energy gap. We find that the emergence of crossing points of ABS is a key 
ingredient to generate sin 2p dependence in the current-phase relation. We further study both the 
energy and p dependence of pair amplitudes in the ferromagnetic region. For large G odd-frequency 
spin-triplet s-wave component is dominant. The magnitude of the odd-frequency pair amplitude is 
enhanced at the energy level of ABS. 


I. INTRODUCTION 

It is known that a number of remarkable physical 
phenomena occur in superconductor/ferromagnet (S/F) 
hybrid structuresFirst one is the generation of tt- 
state^i^ in S/F/S junctions. Since the exchange cou¬ 
pling and spin-singlet Cooper pair are competing each 
other, spin-singlet pairs in ferromagnet have a spatial 
oscillation with changing sign in the presence of the ex¬ 
change coupling A”— Next one is the dominant second- 
harmonic in the current-phase relation of Josephson cur¬ 
rent, sm2p, where p is the phase difference across the 
junction. It is known that sm2p dependenceii appears 
near the O-tt transition point j^^4^ The third one is the 
generation of the odd-frequency pairing in the F region 
by proximity effect in S/F hybrid systemsThere 
have been many theoreticaU^iii and experimentaU^“— 
works about proximity effect via odd-frequency pairing 
in S/F junctions. The fourth one is the so called inverse 
proximity effect where magnetization penetrates into a 
superconductor The electronic property and pairing 
symmetry near the S/F interface is drastically changed 
by this effect. 

Independently of research directions mentioned above, 
study of nanowire on the surface of superconductor in the 
presence of applied Zeeman magnetic field has recently 
become a hot topic in condensed matter physics 
Due to the strong spin-orbit coupling (SOC) in nanowire, 
topological superconducting state is generated. Then, by 
the bulk-edge correspondence, superconducting nanowire 
hosts Majorana fermion (MF) as the end stated”— 
which is one of important factors to realize quantum 


computation4i“— It has also been reported that a chain 
of ferromagnetic atoms on a superconductor forms a 
topologically non-trivial state, where the ABS within the 
superconducting gap is localized around the edge as a 
MFi^— The common feature in these one-dimensional 
topological superconducting systems is that both pair 
potential and magnetization coexist in all sites of the 
nanowire. 

Up to now, although inverse proximity effect and 
topological superconductivity have been studied indepen¬ 
dently, they have not been studied simultaneously. It is a 
challenging issue to clarify a new effect where both effects 
coexist in the same model. If we consider proximity cou¬ 
pled nanowire on the S/F/S junction, we can divide the 
nanowire into the three segments; left superconductor, 
middle ferromagnet and right superconductor. Thus, it is 
possible to design effective one-dimensional S/F/S junc¬ 
tions in nanowire. We consider the situation where the 
ferromagnetic order of F depends on the position in the 
S/F/S junction. Besides, to discuss topological super¬ 
conductivity, we consider Rashba-type SOC in nanowire. 
If ^ which represents the penetration length of ferromag¬ 
netic order is long, zero energy state is generated at the 
S/F interface. On the other hand, if ^ is short, we can 
expect that this zero energy state splits into twc4i^— 
around the magnetic impurity. Thus, proximity coupled 
nanowire on the S/F/S junction is interesting since we 
can study both inverse proximity effect and topological 
superconductivity in the same model. If we tune ^ as 
a parameter, the present model has a unique feature to 
study a new-type of inverse proximity effect including 
topological superconductivity. 
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In this paper, we study electronic spectra and resulting 
Josephson current in this nanowire S/F/S model. First, 
we calculate local density of states (LDOS) of isolated 
left side superconducting segment where magnetization 
penetrates from the right edge proportional to exp[|^^] 
with site index n < 0 and the decay length We clarify 
that if exceeds a certain value, LDOS has a clear zero 
energy peak (ZEP) due to the zero mode at the edge. 
On the other hand, when becomes shorter, LDOS has 
a peak splitting. We call this effect as “zero mode-non¬ 
zero mode crossover”. Throughout this paper, we intro¬ 
duce zero mode and non-zero mode to distinguish zero 
energy state localized at the interface and splitted state. 
We study the Josephson current in the S/F/S nanowire 
junction by changing the decay length of ferromagnetic 
order into the left (right) side superconductor, (Ci?): 
and chemical potential of ferromagnet. It is seen that 
the behavior of Josephson current is quite different when 
the decay length becomes shorter due to the “zero mode- 
non-zero mode crossover”. 

Especially when this non-zero modes^"^— are local¬ 
ized at the boundary between superconductor and fer¬ 
romagnet, we find an anomalous current-phase relation 
which can be roughly expressed as sin2(^. In order to 
understand the physical origin of the current-phase re¬ 
lation more clearly, we calculate the ABS and provide 
an argument that the emergence of crossing points of 
ABS is a key ingredient to produce current-phase relation 
sin 2ip. We find that the emergence of crossing points of 
ABS is a key ingredient to generate sin 2(^9 dependence in 
current-phase relation. This mechanism is distinct from 
preexisting case in S/E/S junctions.— We also calcu¬ 
late pair amplitude decomposing into odd-frequency part 
and even-frequency one. We focus on 5-wave component 
of odd-frequency spin-triplet even-parity (OTE) pairing 
and show their dependence on energy and phase differ¬ 
ence ip. Eor large magnitude of equal spin OTE 

pair amplitude becomes dominant and enhanced at the 
energy level of the ABS. 

This paper is organized as follows: In Sec. im we pro¬ 
vide the formulation to calculate LDOS and Josephson 
current using recursive Green’s function technique. In 
Sec. HITAl we calculate LDOS on the edge of isolated 
left side superconducting segment with decaying ferro¬ 
magnetic order from the right edge. In Sec. ElBl we 
analyze the Josephson current in S/E/S junction chang¬ 
ing decay length of ferromagnetic order parameter and 
chemical potential of ferromagnet. In Sec. IIII Cl we cal¬ 
culate ABS and discuss the relevance to current-phase 
relation in Sec. inrm in sec. itttdi we calculate pair 
amplitudes to shed light on these results from different 
angles. Einally, we summarize our results in Sec. [TVl 


II. FORMULATION 


function. The method of calculating ABS is also included 
in this section. 

First, we review some general aspects of the recursive 
Green’s function technique. Its spirit is as follows: to 
build the full Green’s function, we start from isolated 
blocks and connect them to other parts of the system 
by stacking sites one by one along a certain direction. 
For example, we consider a one-dimensional atomic chain 
along x-direction as shown in Fig. [TJa) and suppose that 
the Green’s function of the detached N sites on the left 
side has been known. We denote ^ as the Green’s 
function at the N site. Gonsider adding one more site 
from the right to this system. Based on Dyson equation, 
we have 


^Ar+i,Ar+i 


■ ^iso T 5^isokw+l, 


( 1 ) 


— ^7V,Ar'^A^W+1^7V+l,Ar+l: 


( 2 ) 


where ^iso is the Green’s function of the isolated N 1 
site. Substituting Eq. ^ into Eq. 0, we get 


^Ar+i,Ar+i 


ftso “ "^AT+l,ATG5(r,AT 

LDOS can be calculated as follows: 


1 -1 


( 3 ) 


p^{E) — ImTr i0~^) 


( 4 ) 


with infinitesimal positive number 0+. Similarly, one can 
determine the Green’s function by stacking from right to 
left (in this case the superscript of G is replaced with 
R, such as G^^). Provided that we start this stacking 
process simultaneously from the two ends, we can obtain 
the Green’s function G^_i at site A" — 1 of the left 
chain and G ^^2 Ar +2 A + 2 of the right chain. 

Erom Eq. (|3]), we know that the process of adding the 


A A + 1 

CK><>0-??-C>0-”0 


(b) N N + l 

CX>11~0""0 ■ 0--OSS-00 

V- 1 N + 2 


FIG. 1: (a) A schematic picture of one-dimensional chain 
which is created by adding sites from left to right, (b) A 
schematic picture of one-dimensional chain which is created 
by stacking sites both from left to right and right to left and 
finally connect the two chains. 

site A in the left chain generates 


In this section we provide a formulation to calculate 
LDOS and Josephson current using recursive Green’s 


G 


L _ 
N,N — 



Vn,n-iG 


L 


Vn-1,N 


-1 


( 5 ) 
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and similarly 

ftso “ ^N-\-1,N-\-2G^_^2,N-\-2^N-\-2,N-\-1 , 

(6) 

from adding the site A/" + 1 in the right chain. Now we 
connect the two chains (light-blue dashed line in Fig. 
[T](b)) and denote the Green’s function of this combined 
chain as G. Based on Dyson equation we get following 
equations: 


/oR _ 

^Ar+i,Ar+i — 


with Boltzman constant temperature T, and Mat- 
subara frequency uOn- 


Next, we construct the model Hamiltonian of the semi¬ 
conducting nanowire on top of the S/F/S junction (Fig. 
[21(a)). We separate this nanowire into three parts by in¬ 
troducing 


Gn,n = — Vn,n-\-iG^_^i^^_^iVn-\-i,n , (7) 

GAr+i,Ar+i = — Vn-\-i,nG^^^Vn,n-\-i (8) 

Gn,n+i = G%^^Vn,n-\-iGn-\-i,n-\-i (9) 

Gn+1,N = ( 10 ) 


51 = {n\l <n< Lsc} 

5 2 = {n\Lsc< n < Lsc + L^m} 

5 3 = {n\Lsc + Lym + 1 < n < 2Lsc + L^m} (12) 


Suppose hopping amplitude of adjacent sites is given by 
t, we can calculate the current 


J = —ieiksTGN-\-l,N{^n) — GN,N-\-l{^n) , (11) 


with site index n. We denote the equal site length of each 
side of superconductors as Lsc and that of ferromagnet 
as I/fm- We define Hamiltonian as follows: 


'H — t cl^^Cno- + ^ ^ + H.C.) flcl^^Cna + H-C-) 


{m,n) 


neSi 

neSs 


neSi 


neSi 


- H MFMcLcna + X] 

neS3 neSs neS2 neS2 


- ^ Vexp 

716*5*1 


n — Lsc 




(clfCnt - cl^Cni) + 


776*53 


Lsc + 7/fm + 1 — n 






(13) 


where c]^^(cncr) is the electron creation (annihilation) op¬ 
erator with site n and spin a, t is the hopping matrix 
between the nearest neighbor (i, j), H is the Rashba spin 
orbit coupling (SOC), /i(/iFM) is chemical potential in 
superconductor (ferromagnet) segment, A is the pair po¬ 
tential, (f is the phase difference of superconductors, W 
is the ferromagnetic order, and ^r(l) is the decay length 
of ferromagnetic order in the right (left) superconductor. 
Unlike the model Hamiltonian on top of the supercon¬ 
ductor with uniform ferromagnetic order, in this model 
construction, we assume that ferromagnetic order pene¬ 
trates into right (left) superconductor segment decaying 
as ^ exp[^](exp[^]) as described in Eq. ([13]). As for 
Rashba SOC, we do not include the second term in Eq. 
m at the interface between superconductor and ferro¬ 
magnet since it does not affect our results. 


Now we apply the recursive Green’s function technique 
to calculate LDOS, Josephson current and ABS. The re¬ 
tarded and Matsubara Green’s function of the isolated 


site can be described as 


9iso{S^) 


1 

{E^io+)-niso 


(14) 


Qisoi^n) — 


1 


-n. 


(15) 


where LLiso is 


'Hiso = ciUi, 


{-jiaoTz - (TyTyA, ne Si 

—n e S2 

-liaoTz - ayTiAe^"^ - ayT^Ae~^'^, n € S3, 

(16) 

with the basis c„ = (c„t, c„^, 4^, 4^). (To,x,v,z{to,x,v,z) 
is Pauli matrix in spin (particle-hole) space and = 
{tx E iTy)l2^ = {tx — iTy)!2. Hopping matrix can be 
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FIG. 2: (a) A schematic picture of a nanowire on top of S/F/S 
junction, (b) Chemical potential of a nanowire in the super¬ 
conductor segment (left) and ferromagnet segment (right). 


written as follows: 






-tcToTz -i^TzCFy, n ^ Si,n ^ S 2 , 
-tcjQTz, n e S 2 


(17) 


III. RESULTS 

This section consists of four parts: in subsection A, 
we study the LDOS on the edge of the nanowire with 
both pair potential and decaying ferromagnetic order and 
exhibit the evolution of the surface resonance modes. 
In subsection B, we then calculate Josephson current 
S/F/S junctions of the nanowire with several different 
decay length. It shows that at specific parameter tuning, 
anomalous current-phase relation which can be roughly 
regarded as sin 2(^9 appears. In subsection C, we calcu¬ 
late ABS. From the spectrum of ABS, we will provide 
a simple argument of explanation that the emergence of 
crossing points of ABS is the important factor to real¬ 
ize the current-phase relation sin2(p. In subsection D, 
we finally analyze the symmetry of pair amplitudes in 
this junction especially focusing on odd-frequency pair¬ 
ing. Odd-frequency spin-triplet pairing is enhanced when 
the Major ana like zero energy state is generated at the 
S/F (F/S) interface's Throughout this section, we fix 
the parameters as follows: A/t = 0.1,///t = = 

'^^Vz/t = 1.5. With this choice of parameters, condition 
of topological non-trivial statoS^ is satisfied when ferro¬ 
magnetic order is uniform. We set the number of sites of 
superconductor segment long enough so that the effect of 
overlapping between two zero energy modes on the edges 
of segment is negligible. In actual numerical calculation, 
we fix the number of sites of superconductor segment as 
Lsc = 4000 and ferromagnet one as Lfm = 4. To cal¬ 
culate retarded Green’s function, we set the infinitesimal 
positive number 0+ as 0+/t = 0.001. 




n-|-l,n — 


-tCFQT- 

-taoTz, n e S 2 


i 2 TzCy^ 


Ti G Si , Ti G S:^ 


(18) 


In the next section, we will calculate the LDOS on the 
edge of the nanowire with proximity coupled pair poten¬ 
tial and decaying ferromagnetic order, Josephson current 
and ABS of the nanowire on S/F/S junction. LDOS is 
given by 


p^{E) = -ImTr 

TT 


Glsc,Lsc(^ + ^0^) 


(19) 


From Eq. Josephson current and ABS are described 
as follows: 


J{ip) = -ieiksT^Tr: G„+i,„(w„)-G'„,„+i(w„) , n e S 2 

( 20 ) 


p{E, ip) =-ImTr Gn,n{E + z0+) 

TT L 


neS2 ( 21 ) 


where (p is the macroscopic phase difference of pair po¬ 
tential between two superconductors. We can calculate 
p dependence explicitly. 


A. LDOS 

In this subsection, we examine the LDOS on the edge 
of nanowire which is calculated from Eq. ([19]) . The situ¬ 
ation is illustrated in Eig. [3]^a), where the pair potential 
and decaying ferromagnetic order coexist. The intensity 
of LDOS is plotted in Eig. [3]^b). It is shown that for the 
long decaying of Vz (^l > 10), there is a single resonance 
peak at zero energy. This is in agreement with the finding 
that in the asymptotic scenario when Vz is spatially uni¬ 
form at infinite such nanowire is topologically non¬ 
trivial which can be confirmed by the sign of Pfafiian with 
our choice of parametersHowever, in the presence of 
spatially nonuniform ferromagnetic order, the preexist¬ 
ing topological argument is no longer valid. Thus, it is 
remarkable to see here that even if the ferromagnetic or¬ 
der is non-uniform, LDOS with zero energy peak can still 
appear on the edge of the nanowire within the numerical 
accuracy. On the opposite, with the decay length less 
than around ^ 10, this zero mode splits into two res¬ 
onance peaks symmetric to zero energy. Owing to the 
localized ferromagnet at the end of nanowire, this find¬ 
ing associates with another important physical results 
known as Shiba statesIt is well known that when 
a magnetic impurity is put on the superconductor, there 
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FIG. 3: (a) A schematic picture of a nanowire on top of super¬ 
conductor with decaying ferromagnetic order, (b) The inten¬ 
sity plot of LDOS on the edge (z.e. n = Lsc) of the nanowire 
shown in (a), (c) The spatial profiles of ferromagnetic or¬ 

der with different decay length, = 25 (black solid line), 
= 5 (red dotted line), = 1 (green dashed line), and 
= 0.2 (blue 2 dotted line), (d)-(g) LDOS on the edge of 
the nanowire with different decay length. 


is a bound state around the impurity inside the super¬ 
conducting energy gap. In the present model, we can 
control the “zero mode-non-zero mode crossover” with 
decreasing the decay length of ferromagnetic order. In 
Figs. [3fd)-(g), we plot the LDOS for a selected decay 
length = 25,5,1, and 0.2 respectively (see also Fig. 
He) where the spatial profiles of ferromagnetic order are 
shown for four different decay length). When = 25, 
we see the zero energy mode, on the other hand, in the 
rest of the cases (Figs. [3fe)-(g)), we find non-zero modes 
in the energy gap of superconducting region. 

In Appendix [Al we consider the overlapping of two zero 
energy modes in the shorter length system and also com¬ 
pare non-uniform ferromagnetic case and uniform case. 


which leads that two zero energy modes appeared in non- 
uniform situation can be regarded as two MFs. 


B. Josephson current 

In this subsection, we study Josephson currents for var¬ 
ious decay lengths of ferromagnetic order on both sides 
of the superconductor. On the left side, we are inter¬ 
ested in three typical decay lengths = 25,1, and 0 
(i.e., ferromagnetic order does not penetrate into the left 
superconductor). In each case, we will tune the decay 
length as well as the chemical potential of the ferro- 
magnet segment {jiiFu/t = —2, —1, 0, and 1, see also Fig. 

Hb)). 


(a) fiFu/t = —2 (b) fiFu/t = —1 



(c) HFu/t = 0 (d) HFu/t = 1 



- a = 25,,eii = 25 

. a = 25,^fl = 5 

___ ^i = 25,£fi = l 
_ a=25,,efl = 0.2 


FIG. 4: (a)-(d) Josephson current of the nanowire on top of 
S/F/S junction in the case of = 25 with four different 
chemical potentials of ferromagnetic layer listed above each 
figure. There are four lines in each figure: — 25, — 25 

(black solid line), = 25, = 5 (red dashed line), = 

25, (fi? = 1 (green dashed line), and — 25, = 25 (blue 2 

dotted line). 

First, we look at the case = 25. As the previ¬ 
ous discussion shows, there is zero energy mode on the 
edge of the left nanowire. In Figs. |4l^a)-(d), we plot 
Josephson currents for /ipivi/^ = —2, —1, 0, and 1, respec¬ 
tively. In each figure, we take four values oi ^r\ 25, 5, 
1 and 0.2. The phase dependence of the current in all 
cases has the dominant coupling proportional to simp. 
Interestingly, for fiFu/t = -2,-1 and 0, the Joseph¬ 
son current in symmetric junctions abruptly changes its 
sign at (/:? = TT as shown in Figs. IH^a), (b) and (c), re¬ 
spectively. In asymmetric junctions, such jump is ab¬ 
sent. Notice that our system considered here has no per- 
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feet transmissivity and is closely similar to that of two 
d-wave superconductor with zero energy ABSs^“— p- 
wave super conductor or Kitaev chains junction sys¬ 

tem where MFs are coupled each otherTherefore, 
the abrupt jump can only be explained by the existence 
of robust zero energy ABSs^“— i.e., MFs^i However, 
in the non-uniform ferromagnetic order, the formation 
of zero modes are distinct from p-wave superconductor 
or Kitaev model. It is interesting that the similar be¬ 
havior of Josephson current is found even in the present 
model. For pFu/t = 1, we obtain 7r-state and the current 
are quite small compared to the other cases of chemical 
potentials. In this case, there is one band at the chemi¬ 
cal potential in ferromagnetic region which has opposite 
spin compared to the ferromagnetic order in SCs. Thus, 
the transparency of the junction is greatly reduced and 
TT-state can appear as a result of misaligned magnetiza¬ 
tions. 


(a) pFu/t = —2 (b) pFu/t = —1 





- - <^L = = 1 

- ^L-h^R-0.2 


FIG. 5: (a)-(d) Josephson current of the nanowire on top 
of S/F/S junction in the case of = 1 with four different 
chemical potentials of ferromagnetic layer listed in above each 
figure. There are four lines in each figure: = 25 

(black solid line), = 5 (red dashed line), = 

l,(fR = 1 (green dashed line), and = 0.2 (blue 2 

dotted line). 

Second, we focus on the case = 1 where non-zero 
modes are localized on the left segment. As for pFu/t = 
—2 and —1 (Figs. [5fa) and (b)), we find the current- 
phase relation sin p in three cases: = 25 (black solid 

line), = 5 (red dotted line), and = 1 (green dashed 
line). On the other hand Josephson current is almost zero 
when ^R = 0.2 (blue 2 dashed line). Surprisingly, when 
PFu/t = —2, ^R = b and = 1 (Fig. [5fc) red dotted 
and green dashed lines), we find the anomalous current- 


phase relation which can be roughly regarded as sin2(p. 
As we will see in the next subsections, the coupling of 
non-zero modes produces this anomalous current-phase 
relation especially when the emergence of crossing points 
of ABS becomes the important factor. For pFu/t = 1, 
all of the Josephson currents are suppressed compared to 
the other cases of chemical potential. 




. ^R = ^ 

= 1 

- ^h = 0.2 


FIG. 6: (a) A schematic picture of a nanowire on top of S/F/S 
junction, (a)-(d) Josephson current of the nanowire on top of 
S/F/S junction in the case of = 0, that is, ferromagnetic 
order does not penetrate into left superconductor segment. 
Results are shown with four different chemical potentials of 
ferromagnetic layer listed above each figure. There are four 
lines in each figure: = 25 (black solid line), = 

= 5 (red dashed line), = 1 (green dashed 

line), and = 0,^i? = 25 (blue 2 dotted line). 

Finally, we study the case = 0, i.e., the ferromag¬ 
netic order does not penetrate into the left supercon¬ 
ductor segment. When pFu/t = ~2, —1 and 1 (Figs. 
[61(b), (c), and (d), respectively), Josephson current is al¬ 
most zero, on the other hand, when fiFu/t = 0, we see 
current-phase relation sin2(^ with decay length on the 
right ^R = 5 (dotted red line) and — sin p with other 
decay length. 

Before we proceed, it is instructive to summarize the 
interesting phenomena we have found in this subsection. 
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In the case of = 25 and fiFu/t = —2, —1, 0, we see the 
abrupt sign reversal of current ai cp = ir. In the case of 
= 25 and pFu/t = 1, the amplitude of Josephson cur¬ 
rent is suppressed. When = 1, P'Fu/t = 0, and = b 
or 1, we obtain current-phase relation approximated as 
sm2(p. We also find this current-phase relation sm2(f in 
the case of = 0, jUFu/t = 0, and ^r = 5. 


C. ABS 

In this subsection, we study the ABS of nanowire on 
S/F/S junction with different decay length and relate 
it to the behavior of Josephson current obtained in the 
previous subsection. We mainly focus on the case of 
I^Fu/t = 0. It is well known that when the magnitudes of 
the pair potential are the same on the left and right side 
of superconductor, Josephson current can be calculated 
by: 


j{^) = 


2e dF{(fi) 
V dip ’ 


( 22 ) 


where F is free energy. In one dimension, F can be writ¬ 
ten as 


F =-fcBry]log(2cosh(^|^)j (23) 

n 

where Sn is the energy of ABS^. Therefore, Josephson 
current is 


J(^) = -^fcBT^ 


cosh(2|fT) 




( 24 ) 


At low temperature, J{(p) can be approximated as 


j(^) -y^sgn(e„) 


dip 


(25) 


sgn(5n) gives -hl(—1) when Sn is positive (negative). In 
the above, n denotes the band index of ABS. Due to the 
particle-hole symmetry, we can only take the ABSs below 
the zero energy into account. Thus, Josephson current 
can be approximated as the derivative of ABSs below the 
zero energy. 

First, we look at the case = ^r = 25 where the zero 
energy modes are located on both sides. These two zero 
modes hybridize as indicated in Fig. [Tl^a). The crossing 
point at (p = TT explains the sudden drop of Josephson 
current which can lead to the unusual Jtt periodicity of 
current-phase relation if we consider AC Josephson cur¬ 
rent. As the decay length on the right decreases, the 
zero energy mode on the left does not hybridize with 
the states on the right, which can be seen from the flat 
ABS as a function of p in Figs. [7fb)-(d). The con¬ 
tribution of Josephson current is mainly carried by the 


PFu/t = 0 

= 25, U = 25 = 25, = 5 



(c) = 25, = 1 (d)a = 25, = 0.2 

1 . 5 . 



FIG. 7: (a)-(d) ABS (the intensity plot of LDOS at the fer- 
romagnet segment within the superconducting energy gap) in 
the case of ppm/t = 0 and = 25 with different decay length 


ABSs away from zero energy. With this estimation, we 
can relate ABS to the behavior of Josephson current 
in the previous subsection. If we look at Figs. m 
and (c) ((d)), ABSs below the zero energy change as 
^ — cos p (^ cos p) which leads to current-phase rela¬ 
tion sinp (—sinp). This consideration corresponds well 
to the black solid line and red dotted line (green dashed 
line) in Fig. SJc). If we tune the chemical potential of 
ferromagnetic layer as pFu/t = 1, the ABSs are almost 
flat which do not make contribution to the Josephson 
current (Fig. [8j). Indeed, the amplitude of Josephson 
currents are almost zero as shown in Fig. IH^d). 

Next, we focus on the case = 1, when non-zero 
modes are localized at the interface between left super¬ 
conductor and ferromagnet segment. When decay length 
^R = 25 (Fig. [9l^a)), the zero energy mode is located on 
the right segment. This mode is not coupled with the 
left non-zero modes. The major change of ABSs below 
the zero energy can be regarded as ^ — cos(/:? (Fig. (Ifa)) 
which contributes smp to the Josephson current (Fig. 
[5](c) black solid line). 

In the case of ^r = 5 and = 1, when the anomalous 
current-phase relation sm2p can be seen, we find that 
non-zero modes on the both sides hybridize and these 
states are crossed at two values of p: one is located at 
between p = 0 and tt, another is between p = tt and 
p = 27r (Figs. Mjc>) and (c)). As we discuss later, this 
crossing points of ABS is important to realize current- 
phase relation sin2(^. For ^r = 0.2, all of the ABSs 
change as cosp (Fig. M^d)) which gives current-phase 
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(a).. ^ 


/^fmA — 1 

= 25 (b) fr =: 


£L = 25,eR = 25 V^;_.eL = 25,efl = 5 



(c) iL = 2&fl = 1 = ‘^i^R = 



/xfmA — 0 

= 1,^/J = 25 = 

l.M 


= 5 




- 1.51 _ ^ 

0 1 i 0 1 

CL = l,^ii = l = = 

1 . 5 | 



FIG. 8: (a)-(d) ABS (the intensity plot of LDOS at the fer- 
romagnet segment within the superconducting energy gap) in 
the case of /xpivi/t = 1 and = 25 with different decay length 
^R- 


FIG. 9: (a)-(d) ABS (the intensity plot of LDOS at the fer- 
romagnet segment within the superconducting energy gap) in 
the case of iiYm/t = 0 and = 1 with different decay length 
^R- 


relation simp (Fig. [5fc) blue 2 dotted line). 

Finally, we focus on the case = 0, Le., ferromagnetic 
order is absent in the left superconductor side (Fig. [TQ|) . 
If we set = 25, all of the ABSs becomes flat as a func¬ 
tion of (p which leads to almost zero Josephson current 
(Fig. [bfc) black solid line). For = 1,0.2, the change 
of ABS below the zero energy obeys ^ cos(p (Figs. flOlc) 
and (d)), thus, current-phase relation is — siiKp (Fig. [6fc) 
green dashed line and blue 2 dotted line). On the other 
hand, in the case of = 5, when we see current-phase 
relation sm2(p, we again find two crossing points of ABS 
(Fig. Mh))- 

From these analysis explained above, we find that the 
behavior of Josephson current can be determined by the 
hybridization of ABSs and when this hybridization pro¬ 
duces the two crossing points at the zero energy, current- 
phase relation sm2(p can be seen. Now we consider the 
simplified situation where there are two crossing points 
of ABSs: we have only two states within the energy gap 
(blue area in Fig. [TT]top) which are labeled by ea and sr 
(Fig. [TT]top). We assume sa (^b) is written as ^ — cos(^ 
(cos (p) and crossing points are located at (p = ^ and ^. 
According to Eq. (|25|) and particle-hole symmetry, we 
can focus on the states below the zero energy to esti¬ 
mate Josephson current. We separate the region of cp 
into three: I {(p\0 < (p < f}, II {^|f < ^ < 

III {(p\^ < ^ < 27r} (Fig. [11] top). In I and III, 
ABS below the zero energy obeys — cos(/9, thus, Joseph¬ 
son current reads ^ sin(^, while in II, ABS transforms 
into simp which gives Josephson current —simp. We 


plot Josephson current as a function of (p {0 < (p < 27r) 
in bottom of Fig. [TTl Due to the different curve that 
ABS obeys in I and II (II and III), there is a jump at 
the boundary between I and II (II and III) which gen¬ 
erates current-phase relation sin2(^. Up to now, sin2(^ 
dependence of Josepshon current has been discussed in 
S/ferromagnetic insulator/S junction^ S/F/S with dif¬ 
fusive F near the vicinity of 0 — tt transition point 
d-wave superconductor junctions^i^“— and s-wave/spin 
triplet p-wave superconductor junctionsOur setup 
of realizing sin 2 (p dependence is distinct from the preex¬ 
isting cases. 


D. Symmetries of Cooper pair 

In this subsection, we focus on the symmetry of 
Cooper pair in the present one-dimensional S/F/S junc¬ 
tions. In odd-frequency pairing state, the pair amplitude 
changes its sign with the exchange of times of two paired 
electronsTaking account of this symmetry class, sym¬ 
metry of Cooper pair is classified into (1) even-frequency 
spin-singlet even-parity (ESE), (2) even-frequency spin- 
triplet odd-parity (ETO), (3) odd-frequency spin-triplet 
even-parity (OTE)^ and (4) odd-frequency spin-singlet 
odd-parity (OSO)^^ Although odd-frequency bulk su¬ 
perconductor has not been discovered up to now, odd- 
frequency pair amplitude can exist ubiquitously as a sub¬ 
dominant state. It is known that odd-frequency pairing is 
induced by the breaking of the translational^”— or spin- 
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fJ'Fu/t = 0 

(a) eL = 0,ei? = 25 (b)a = 0,^fl = 5 



FIG. 10: (a)-(d) ABS (the intensity plot of LDOS at the fer- 
romagnet segment within the superconducting energy gap) in 
the case of /j^Fm/t = 0 and = 0 {i.e. there is no ferromag¬ 
netic order in the left superconductor segment.) with different 
decay length ^r. 


rotational symmetry^ i^^i^^ from bulk even-frequency pair 
potential. Also, it has been clarified that zero-energy lo¬ 
cal density of states are enhanced by the odd-frequency 
pairingOdd-frequency pairing influences seri¬ 
ously the proximity effect and various electronic proper¬ 
ties of the junctions. 

In the present nanowire S/F/S junctions, the symme¬ 
try of pair amplitude far from the S/F (F/S) interface 
is ESE s-wave one, since the induced pair potential is 
conventional spin-singlet 5-wave. Due to the symmetry 
breaking, near the S/E (E/S) interface or inside ferro- 
magnet region, odd-frequency pairings can be induced. 
Here, we focus on 5-wave component of ESE and OTE 
pair amplitudes in ferromagnet region. Eirst we focus 
on the real frequency representation of pair amplitudes. 
ESE and OTE amplitudes are given by 

feveniE) = ^+ i0+) + F^i{-E - i0+)} (26) 



FIG. 11: (top) A schematic picture of two ABSs denoted by as 
(red) and sb (green). These two states have two crossing 
points (solid circle) at E = 0. (bottom) Josephson current 
which is given by J((p) ~ — ■^^sgn(£A) 


that there is a close relation between ABS and odd- 
frequency pairing. In the presence of zero energy ABS as 
a surface state of unconventional superconductors, odd- 
frequency pair amplitude is hugely enhanced i^^i^^F^T^ 
Thus the presence of zero energy state (ZES) can be inter¬ 
preted as an emergence of odd-frequency pairing. Also, 
it has been clarified that ME always accompanies odd- 
frequency pairing 


a = 25,^i?, = 25, /xfm/^ = 0 

(a) Re[OTE;|] (b) Im[OTE|;] 



foM (E) = \ {E + i0+) - Elii-E - *0+)} (27) 
with n G 5'2- comes from Eq. ©: 


G 


n,n — 


G E 
F G 


( Ft.t pt.-l- 

r _ / n,n ^ n,n 

~ I z^t,t 

\ n,n n,n 


(28) 


Above representation is useful to compare with the ABS. 
To clarify the relation between the ABS and the pair am¬ 
plitude is an interesting issue, since it has been revealed 


FIG. 12: (a)(b):The intensity plot of real part (a) and imagi¬ 
nary part (b) of odd-frequency spin-triplet even-parity (OTE) 
pair amplitude with H spin component fodtiiF) in the case 
of = 25, = 25, fiFm/t = 0. 


Eirst, we calculate pair amplitude for = 25, = 25 

and = 0. The corresponding ABS with the hy¬ 

bridization of Majorana like ZESs in left and right side 
has been shown in Eig. [Tl^a). In this case, OTE pairing 
becomes dominant for any phase difference. As compared 
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to other components, feven{E), the 

magnitude of fodd^E) is dominant. This is because that 
the direction of the majority spin in ferromagnet is down. 
Therefore, we plot s-wave OTE pairing f^dd^E) in the 
middle of the ferromagnet in Fig. [121 The energy and 
Lp dependence of is almost similar to those of 

ABS in Fig. [71 For (/:? = 0, tt and 27r, Re[f^^^{E)] = 

-Mfhdi-E)] andIm[/,Y,(E)] =Im[/„Vd(-^)]- These 
relations of odd-frequency pairing are known in the previ¬ 
ous study in normal metal / superconductor junctions 
Next, we look at the case with ^ and fipM = 

0. The corresponding ABS with double crossing points 
in the energy spectrum of ABS has been shown in Fig. 
He). Such characteristic also appears in the intensity 
plot of ftaiiE), PJaiE) and feven{E). The E and cp de¬ 
pendences of fo^^{E) are similar to those of fo^{E). The 
remarkable point here is that not only odd-frequency pair 
amplitude but also even-frequency pair amplitude exists 
with the same order in contrast to the case in Fig. [131 
For ^ = 0, TT and 27r, ReiOUE)] = -Re[f^fP-E)] and 
(^)] = (-^)] with cr =t 4) and a' =t 

(4,). On the other hand, Re[/e«en(-E’)] = Re[/e«en{--E’)] 
and Im[/ocid(R')] = -'i-^[fodd{-E)] are satisfied. 


IV. SUMMARY AND CONCLUSIONS 

In this paper, we have studied LDOS, current phase 
relation of Josephson current, energy levels of Andreev 
bound state and induced odd-frequency pairings in su¬ 
perconductor / ferromagnet / superconductor nanowire 
junction, where the magnetization penetrates into super¬ 
conducting segment with a decay length We have cho¬ 
sen the chemical potential and SOC so that the topologi¬ 
cal superconducting regime hosting MF is realized for suf¬ 
ficiently large magnitude of We have found that when 
^ becomes larger, LDOS has a ZEP. On the other hand, 
if ^ is shorter, zero energy state at the interface between 
superconductor and ferromagnet splits into two states. 
Accordingly, the behavior of Josephson current drasti¬ 
cally changes. By tuning the parameters of the model, 
we have found an almost second-harmonic current-phase 
relation, sin2(/9, with phase difference cp. Based on the 
analysis of ABS, we clarify that current-phase relation is 
determined by coupling of the states within the energy 
gap. We find that the emergence of crossing points of 
ABS is a key ingredient to generate sin2(^ dependence in 
current-phase relation. We further studied both the en¬ 
ergy and cp dependence of pair amplitudes in ferromagnet 
region. For long odd-frequency 5-wave triplet compo¬ 
nent is dominant. The magnitude of the odd-frequency 
pair amplitude is enhanced at the energy level of ABS. On 
the other hand, when ^ becomes shorter, not only odd- 
frequency pairing but also even-frequency pairing mixes. 

Recently, sin2(/? behavior has been observed in S/F/S 
junctions^ when ferromagnet is an insulator which has 


= 1 5 f^Fu/t = 0 

(a) Re[OTEn] (b) Im[OTEtJ 



(c) RelOTE^J (d) Im[OTE|^] 



FIG. 13: (a)-(d): The intensity plot of real part (a) ((c)) and 
imaginary part (b)((d)) of odd-frequency spin-triplet even- 
parity (OTE) pair amplitude with spin component ftkiiE) 
and H spin component foM(E)- In (e) and (f), real part and 
imaginary part of even-frequency spin-singlet pair amplitude 

feven{E) IS plotted. 


a spin filter effect. Thus to clarify the relevance of our 
obtained sin2(/:? dependence to this experimental report 
is a challenging issue. In our paper, ballistic transport 
is assumed. If ferromagnet becomes diffusive, we can 
expect anomalous proximity effec t by odd- 
frequency spin-triplet 5-wave pairing. Extension to this 
direction is also an interesting future study. 
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Appendix A: Spatial profile of LDOS 

In this appendix, we show spatial profile of LDOS 
of superconductor with decaying ferromagnetic order to 
understand “zero mode-non-zero mode crossover” more 
clearly. The ferromagnetic order decays from right to left 
similarly to Eig. Ha). Since we will consider hybridiza¬ 
tion of two ZESs on the both edges of superconductor, 
we set the number of sites rather short, 200. 

These profiles are shown in Eig. [TH Eor long enough 
decay length (Eig. HH^a)), ZEPs are localized on the both 
side of the nanowire, which is analogous to the situation 
where MEs are localized on the both edges of topolog¬ 
ical superconductorIf the decay length is decreased, 
however, the zero energy mode on the left comes close to 
the that on the right edge (Eig. flTlbB. then finally these 
two zero modes hybridize (Eig. [Mf cB to split into two 

(Fig. Hd)). 


(a) a = 625 



(b) a = 125 




FIG. 14: (a)-(d) The intensity plot of LDOS of the nanowire 
with decaying ferromagnetic order on top of the supercon¬ 
ductor. The ferromagnetic order decays from right to left 
similarly to Fig. [3l^a). Horizontal axis represents site index 
and vertical one does energy. Decay length of the ferromag¬ 
netic order is set as (a)^L = 625, (b)^L = 125, (c)^l = 25, 
and (d)^L = 0.2. 



Noriiiiil Nanowire+SC 
metal +ciec(n ing FM 


FIG. 15: (a)-(c) The spatial dependence of LDOS at zero en¬ 
ergy of the semiconducting nanowire with pair potential and 
decaying ferromagnetic order whose length is 200 sites. The 
several cases of LDOS are plotted for different decay length 
which is indicated by different colors as shown in each fig¬ 
ure. (d) The spatial dependence of LDOS of the normal 
metal/nanowire junction system. The total length of this 
junction is set as 200 sites and the interface between normal 
metal and nanowire is located at site 140. In the nanowire 
segment, the pair potential and decaying ferromagnetic or¬ 
der with = 125 are also included similarly to (a)-(c). A 
schematic picture of this junction is shown below the plot. 


Below, we analyze the spatial profile of LDOS of the 
system we consider here in more detail especially focus¬ 
ing on at the zero energy. In Figs. [TSl aVfcL we plot this 
spatial dependence of LDOS for several cases of decay 
length. For large enough (Fig. [T5l aj and (b)), we see 
the two ZEP: one is on the edge, i.e. at 200 sites and an¬ 
other which is identified as broad peak is left from it (for 
Fig. [T5lf aL this ZEP is located at site 140 ^ 150 and for 
Fig. [T5lbL it is at around site 180). We notice that the 
position of the left peak is shifted if the decay length is 
changed and accordingly the behavior of the oscillations 
of LDOS left from the broad peak changes. However, we 
also find that the oscillations between the two peaks show 
the similar behavior. When the decay length is shorter 
than around = 20, the two peaks hybridize and are 
away from the zero energy (Fig. fTSf cB. Eurther, we 
also calculate the spatial dependence of LDOS of normal 
metal/nanowire junction system with = 125. We set 
the interface at the site 140. The result is shown in Eig. 
[T5lf dL The board peak which is indicated in Eig. [T5f ai 
spreads into the normal metal. 
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(a) 



Normal Nanowire+SC 
metal ^uniform FM 


(b) 

12 



Nanowire+SC 
-^-uniform FM 


FIG. 16: (a) The spatial dependence of LDOS of the normal 
metal/nanowire with pair potential and uniform ferromag¬ 
netic order junction system whose length is 200 sites. The 
interface is positioned at 50 sites, (b) The similar plot to (a) 
but the normal metal is replaced by s-wave superconductor. 
Schematic pictures of the junction are shown below each plot. 


calculate the spatial profile of the LDOS of the normal 
metal/nanowire with pair potential and uniform ferro¬ 
magnetic order junction system as well as that of s-wave 
superconductor/nanowire junction system. These plots 
are shown in Figs. M. a) and (b), respectively. The in¬ 
terface is located at site 50. In Fig. [T5la). the ZEP on the 
left which should be expected to appear at the interface 
penetrates into the normal metal segment. This is anal¬ 
ogous to Fig. fl^ d). Moreover, in Fig. flGr b). we can see 
the ZEP on the left (at site 50), however, as opposed to 
Figs. flCTa). this ZEP cannot spread into the left from the 
interface due to the superconducting energy gap. There¬ 
fore, ZEP cannot be stabilized and decays away from the 
interface. This behavior is similar to Eigs. [T^a)(b). how¬ 
ever, there is one difference: in Eigs. [T5fa)(b) due to the 
presence of non-uniform ferromagnetic order, there isn’t 
the explicit boundary which distinguishes the topologi¬ 
cally trivial area and non-trivial one. Thus, the broad 
ZEP appears in Eigs. [T^a)(b). 


Actually, we can regard the two ZEPs explained above 
as the MEs. To make this statement more convincing, we 


The comparison between the results of Eig. [15] and 
Eig. [T6l implies that two ZEPs appeared in non-uniform 
case can be identified as two MEs. 
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